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ABSTRACT 



We examine the magnetic field in NGC 4214, a nearby irregular galaxy, using 
multi-wavelength radio continuum polarization data from the Very Large Array. 
We find that the global radio continuum spectrum shows signs that free-free ab- 
sorption and/or synchrotron losses may be important. The 3 cm radio continuum 
morphology is similar to that of the Ha while the 20 cm emission is more dif- 
fuse. We estimate that 50% of the radio continuum emission in the center of the 
galaxy is thermal. Our estimate of the magnetic field strength is 30 ± 9.5 fiG 
in the center and 10 ± 3 fiG at the edges. We find that the hot gas, magnetic, 
and the gravitational pressures are all the same order of magnitude. Inside the 
central star forming regions, we find that the thermal and turbulent pressures 
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of the HII regions dominate the pressure balance. We do not detect any signifi- 
cant polarization on size scales greater than 200 pc. We place an upper limit of 
8 /j,G on the uniform field strength in this galaxy. We suggest that the diffuse 
synchrotron region, seen to the north of the main body of emission at 20 cm, 
is elongated due to a uniform magnetic field with a maximum field strength of 
7.6 /LtG. We find that, while the shear in NGC 4214 is comparable to that of 
the Milky Way, the supernova rate is half that of the Milky Way and suggest 
that the star formation episode in NGC 4214 needs additional time to build up 
enough turbulence to drive an a — u dynamo. 

Subject headings: galaxies: ISM - galaxies: individual (NGC 4214) - galaxies: 
irregular - galaxies: magnetic fields - polarization 



1. Introduction 



Modern theories of galaxy formation postulate a bottom-up process where smaller sys- 
tems merge together to for m successively larger structures like the spiral galaxies we see 
today (jSpringel et al.ll2005l ). In the local universe, low-mass galaxies are analogs to these 
high-redshift "building blocks." By studying nearby low mass galaxies, in particular low 
mass irregular galaxies, we can gain detailed insight into the basic physical processes at 
work in these systems, which are different than those in spirals or ellipticals due to their 
smaller mass. Studying the interstellar medium (ISM) of irregular galaxies is particularly 
important because their ISM is strongly influenced by star formation and interactions, which 
significantly affect the future evolution of these galaxies. 

There have been a number of surveys of low mass irregular galaxies aimed at investigat- 
ing the properties of the neutral and ionized g as components of the ISM and the contributions 
to the ISM by the stellar components (e . g.. iHunter fe Elmegreenl I2004J : iHunter et al.l 12006 



Hunter fc Elmegreenl l2006t IHunter et al.l 120071 . 120101 ) . Relatively little attention, however, 
has been paid to the magnetic field component of the interstellar medium. 

In general, the magnetic field energy is thought to be in equipartition wi th the turbu- 
lent and cosmic ray ene rgies (jZweibel &: Heilesl 119971 ; IFerrierd l2001r I Cox! l2005t although see 
Chi fc Wolfendald Il993l for a contrasting view) and is thus an important par t of the inter- 
stellar med i um. In the few irr e gular galaxies with measured magnetic fields (IHavnes et al. 



1986 



199 



Mao et al. 



Chyzv et al. 



2008 



2000 



Chyzv et al. 



2000: 



Chyzv et al. 2003: Gaensler et al. 



20051 : 



Kepley et al.l I201CH ; IChyzy et al.l 120111 ) , th e magnetic fie l ds are comparable 



in strength to those in larger spirals like the Milky Way ( IFerrierd l200ll : IWielebinskil 120051 ) 
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and M51 ( IBeck et al.lll987l ). Therefore, in irregular galaxies, the magnetic field may play a 
more important role in the dynamics of the ISM than in spir al galaxies. Magnetic fields also 



chann el gas flows and distribute and accelerate cosmic rays (IBeckl l2005r Everett et al.l 12008 



2010l ) . Irregular galaxies have shallow potential wells that are easily dis rupted and poten- 
tially have a significant am ount of their mass expelled from the galaxy (IFerrara fc Tolstoy 
2 OOOl ; iTremonti et al.l 12004 ) . Magnetic fields may either help the baryons e scape or help to 



confine them to the galaxy depending on the structure of the magnetic fields (lTomisakalll990 



Mineshige et al.lll993 : Ide Avillez & Breitschwerdtl 120051 : Everett et al.ll2008l . l2010f ). Finally, 
the presence of significant magnetic fields in irregular galaxies raises the question of how they 
were generated in the first place. The main ingredients in generating and sustaining galactic 
magnetic fields are generally thought to be large-scale shear and small-scal e turbulence. The 
well known a — uj d ynam o models are based on the se effects (although see IZweibel fc Heiles 



19971 . lKulsrudlll999l . , and Kulsrud fc Zweibelll2008l for critiques of this model). In the a — u 
dynamo, differential rotation (cj) is used to stretch small-scale turbulence (a) leading to an 
increase in magnetic field strength and the presence of a large-scale magnetic field. It may 
be harder for a — u dynamos to operate in irregular galaxies because they rotate mostly as 
solid bodies and their ISM is more easily disrupted because of their smaller potential wells. 
The strength and structure of the magnetic fields in irregular galaxies provide stringent 
constraints for dynamo models of these objects. 

At pre sent, only a handful of nearby galaxies have detailed magnetic field measurements: 



1991 



NGC 1569 (IKeplev et al. 



Klein et al 



1993 



2010). the Large and Small Magellanic Clouds (IHavnes et al 



Gacns 



er et al. 2005 



2000), IC 10 flChvzv et al.l 120031 ). NGC 6822 flChvzv et al.ll2003l ). and IC 1613 flChvzv et al. 



Mao et al.l 12008). NGC 4449 fC hvzv et al. 



1986 



201 ll ). To better understand the role of magnetic fields in the interstellar medium of irregular 
galaxies and how these fields are generated in these systems, we have embarked on a program 
to obtain deep, high-res olution radio contin uum data on several irregular galaxies. In this 
first paper in this series (IKepley et al.ll2010l ). we presented our observations of the magnetic 
field of NGC 1569, a dwarf irregular galaxy hosting one of the most extreme starbursts in 
the nearby universe. The magnetic field in this galaxy is shaped by its outflow and the field 
may be playing an important role in channeling gas away from the disk. 

Work by other groups have found that the LMC, SMC, and NGC 4449 have large-scale 
ordered fields like NGC 1569. The LMC has a 1 11 G regular, axisymmet ric spiral field, a 
random field of 4.1 /iG, and a total field of 4.3 /xG ( jGaensler et al.l 120051 ) . A cosmic ray- 
driven dynamo ( )Parkerlll992t lHanasz et al.l 120041 ) has been invoked to expla in the presence 



of a regular field despite the vigorous recent star formation in the LMC (jGaensler et al. 



20051 ). The SMC has a random magnetic field strength of 3 fiG and a unidirectional, regular 
field with a strength of 1.7 fiG aligned roughly along the Magellanic bridge, which suggests 
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an origin associated with the interaction between the LMC and SMC (IMao et al.l 120081 ). 
In contrast to the weak, but regular fields in the L MC and SMC, NG C 4449 has a strong 
uniform field of 8 fiG and a total field of 14 /zG ( IChyzy et al.l 120001 ) . The uniform field 
shows a fan-like structure in the inner regions of the galaxy and a spiral structure in the 
outer regions. A fast dynamo model including velocity shear around a bar, a wind from the 
central region of the galaxy, a nd a random field produces a ma gnetic field configuration very 
similar to that of NGC 4449 (lOtmianowska-Mazur et al.ll2000l ). In contrast to the magnetic 
fields of NGC 1569, the LMC, the SMC, and NCG 4449, the uniform fields of IC 10, NGC 
6822, and IC 1613 (as trace d by the polariza tion of their radio continuum emission) are 
either weak or non-existent ( IChyzy et al.l 120031 ) . IC 10 has a diffuse, random field with a 
strength of 14 /j,G in the center and 7 /zG on its outskirts. NGC 6822 has a random field 
with a maximum strength of 5 /iG and IC 1613 has a random magnetic field of 2.8 /j,G. Both 
IC 10 and NGC 6822 have small patches of unif orm fields with stren gths of 2-3 /zG on size 



scales much smaller than the size of the galaxy. IChyzy et al. (2003) sugges t that the fields 
in these galaxies are the result of a fluctuating dynamo (ISubramanianlll998l ). IC 1613 shows 
no sign of uniform field; any polarized emission can be attributed to background sources 
rtchvzv et al-lboilh . 



The line of sight fields determined by lGaensler et al.l ( 120051 ) and lMao et al.l (120081 ) for the 
LMC and SMC are based on the rotation measures of background sources, while the results 
for NGC 4449, IC 10, NGC 6822, IC 1613, and NGC 1569 are based on observations of the 
radio continuum polarization. There are significant differences between the two methods. 
The radio continuum polarization observations trace the magnetic field with cosmic ray 
electrons, while the rotation measure observations trace the field with thermal electrons. In 
addition, only the rotation measure observations can tell a regular (i.e., coherent) field (all 
one direction) from a uniform field (one orientation, but the direction uncertain by 180°). 
Unfortunately, current technology limits this method to nearby galaxies (e.g., LMC, SMC, 
M31) with very large angular siz es, although this situati on will improve with the development 
of the Square Kilometer Array ( IBeck fe Gaenslerll2004j ). 



In this paper, we foc us on NGC 4214, whi c h is c lassified as a mixed, S-shaped, Mag- 



ellanic irregular galaxy by Ide Vaucouleurs et al.l (119951 ). NGC 4214 is a well-studied galaxy 
whose rest-frame UV spectrum strongly resembles that of star-forming galaxies at z = 3 



(Stcidcl et al. 


1996). 


Observations of its stars and ionized gas (' 


Vtaiz-Apellaniz et al. 


1998; 


MacKentv et al. 


2000 


Drozdovskv et al. 


2002; 


Ubeda et al. 


2007a 


b ) show that there are two 



active regions of star formation in the galaxy: NGC 4214-1 (also referred to as the northwest 
complex) and NGC 4214-11 (also referred to as the southeast complex). NGC 4214-1 is the 
older co mplex with an age of 3 to 4 Myr and consists of an Ha shell surrounding a super-star 



cluster (IMacKenty et al.ll2000f ). The Ha kinematics of NGC 4214-1 show some indications of 
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expanding shells (IMartmlll998l : iMafz-Apellaniz et al.lll99a IWilcots fc Thurowl 120011) . NGC 
4214-11 is slightly younger with an age of 2.5 to 3.0 Myr (IMacKenty et al.ll2000f ). The Ha in 
this region is coincident with the sta r cluster, showing that t he cluster has not had suffic ient 
time to clear out a large Ha cavity (IMacKenty et al.ll2000i ). ISchwartz & Martini (120041 ) de- 
tect an outflow of cold gas in this region with a velocity of 23 km s -1 , which is associated 
with a small Ha bubble. X-ray observations with Chandra an d XMM- Newton show that the 



hot ga s in this galaxy is centered on the star forming regions (IHartwell et al.ll2004t lOtt et al. 



2005cM . The star forming regions of NG C 4214 and their a s sociated ionized and hot gas are 
embedded in a disk of much older stars (jFanelli et al .1119971 : iDrozdovsky et al.l 120021 ). which 
in turn is embedded in a neu tral hydrogen e nvelope which e xtends to 1.4 times the Holmberg 
radius (5.3'fj of this galaxy jAllsopplll979l : lMcIntvre!ll997h . 



The dist ance to NGC 4214 is controve rsial, with dist ances used in the literature ranging 
from 2 Mpc JHopp et al.lll999r i to 7 Mpc (Ullsoppl \l97$ \ . We adopt a distance of 2.94 Mpc 
in this paper, which is the distance for the galaxy d etermined from Hubb l e Spa ce Tele- 
scope imaging of the tip of the red giant branch by IMafz-Apellaniz et al.l ( 120021 ). This 
dis tance agrees with the ti p of the red gi ant branch distances d etermined independently 
bv brozdovskv et al.l (hood . 2.70 Mpc) and balcanton et all (hood . 3.04 Mpc) and w i th the 



distance determined using the planetary nebula luminosity function by iDopita et al.l (12010 



3.19 Mpc). It is also in reasonable agreement with distances determined by its radial velocity 
and models of the Hubble flow determined by NEE0 (3.75 to 7.53 Mpc). At the adopted 
distance, 1" is 14.2 pc. 

The observations presented in this paper represent the most sensitive observations of 
the radio continuum emission of NGC 4214 to date. Previous observations of the radio 



contin uum of NGC 4214 include lAllsoppI (119791 ) . iBeck et all (120001 ). an d IMacKenty et al. 
( 2000). Although the resol ution of our observations is lower than that of Beck et al. ( 200oh 



and lMacKenty et al.l (120001 ). we integrated for substantially longer (effective integration times 
ranging from 5 hours to 15 hours rather than 15 minutes). Our choice of array configuration 
(see Section [2]) allows us good sensitivity to both small scale features and large-scale diffuse 
emission in NGC 4214. Our observations trace the thermal emission produced by the ionizing 
radiation of young stars and the synchrotron emission produced by the detonation of old 
massive stars, which nicely complements the new data on the stellar content of NGC 4214 



obtained by the Wide Field Camera 3 on the Hubble Space Telescope (e.g., IDopita et al. 



1 The Holmberg radius is the radius where the surface brightness of a galaxy falls to 26.5 mag arcsec 2 . 

2 Thc NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology, under contract with the National Aeronautics and Space Administration. 
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20101 ) and the new data o n the molecular components of the ISM obtained using Herschel 



.e.g., 



Cormier et al.ll2010f ). 



We detail our data calibration and imaging process in Section [2J In Section [31 we 
present our observations in detail. Section 0] is devoted to our results including a comparison 
with multi- wavelength data (Section 14. ip . estimates of the thermal fraction of the radio 
continuum emission (Section 14.21) . an estimation of the magnetic field strength in NGC 4214 
(Section I4.3p . an analysis of the importance of the magnetic field in the ISM of NGC 4214 
(Section l4.4p . an investigation of the magnetic field structure of NGC 4214 (Section 14. 5p . 
estimates of cosmic ray lifetimes in NGC 4214 (Section 14 . 6 j) . and a discussion of possible 
mechanisms for generating the field (Section 14. 7p . Finally, we summarize our conclusions in 
Section [5J 



2. Data 

We have obtained radio continuum polarization observations of NGC 4214 at 20, 6, and 
3 cm using the National Radio Astronomy Observatory (NRAO) Very Large Array (VLA)E 
Our observations (proposal IDs: AK606, AK616; PI: Kepley) were taken between 2005 
August and 2005 December. See Tabled] for a summary of these observations. We used two 
different configurations of the VLA to as closely as possibly match the synthesized beams at 
our three different observing frequencies. The correlator was configured for two intermediate 
frequencies (IFs) each with 50 MHz bandwidth and full polarization products (RR, LL, 
RL, LR). For all our observations, we used 1331+305 (3C286) to calibrate the flux density 
scale and absolute polarization angle and the relevant secondary calibrator to calibrate the 
amplitudes and phases as well as the instrumental polarization. 

For our 20 cm observations, we observed our secondary calibrator (1227+365) for four 
minutes for every thirty-five minutes on source. The primary calibrator (1331+305) was 
observed for ten minutes near the beginning of the observations and ten minutes at the end 
of the observations. For our 6 cm observations, we observed each pointing of a four pointing 
mosaic for eight minutes and observed the secondary calibrator for four minutes once every 
cycle through the mosaic. The primary calibrator was observed in a fashion similar to 
that of our 20 cm observations. For our 3 cm observations, we observed each pointing of a 
sixteen pointing mosaic for six minutes. Every half cycle through the mosaic we observed 
our secondary calibrator (1146+399) for six minutes. The primary calibrator observations 



3 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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were the same as in our 20 cm and 6 cm observations. 

We calibrated the data in AIPS. Our data reduction procedu res followed th e standard 



procedures to calibrate VLA data outlined in the AIPS Cookbook (lGreisenll2008l ) with some 
additional effort to properly calibrate the polarization. We briefly sketch our data reduction 
procedures below. 

After the initial flagging of the data, we used the primary calibrator (with the appro- 
priate u-v limits) to set the primary flux density scale. With our observational setup, 3C286 
(1331+305) is well-approximated as a point source. We then calibrated the amplitudes and 
phases of the primary and secondary calibrators. The flux density value for the secondary 
calibrator was bootstrapped from the flux density value obtained for the primary calibra- 
tor. We applied the calibrated amplitudes and phases to their respective calibrators and 
the calibrated amplitudes and phases from the secondary calibrator to the source. To cal- 
ibrate the instrumental polarization response, we used our observations of the secondary 
calibrator. The secondary calibrator was observed at least five times over at least 100 de- 
grees of parallactic angle, which allows us to use the rotation of the feed with respect to 
the sky to separate the instrumental polarization (which does not rotate as the parallactic 
angle changes) from the astronomical polarization (which does rotate as the parallactic angle 
changes). The calibration for the instrumental polarization was then applied to the data. 
The primary calibrator was used to determine the absolute polarization angle and this final 
correction applied to the data. 

After calib rating in AIPS, we exported the I, Q, and U data as FITS files and imaged the 



data in Miriad (ISault et al.lll995l ) to take advantage its ability to jointly deconvolve mosaics. 
For the 3 cm and 6 cm data, we used the mosaic option on the invert task to combine the 
different pointings of the mosiac. The data was weighted with a robust value of 0. For 
the 3 cm data, we tapered the data in the u-v plane to match the resolution of the data at 
20 cm and 6 cm. We then cleaned the mosaic using the task mossdi. Since the 20 cm data 
was only a single pointing, we used invert without the mosaic option to image the data and 
the task clean to clean the resulting image. In the 20 cm data, there was a bright source 
to the southwest of the galaxy that we were not able to completely model during cleaning. 
The incomplete removal of the source has the side-effect of increasing the noise in the 20 cm 
Stokes I image to approximately 4 times the theoretical noise at 20 cm and twice the noise 
in the other images. The Stokes Q and U images were not affected as severely and their 
noise levels are much closer to the theoretical noise. Table |2] summarizes the properties of 
our final images. 
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3. Properties of the Radio Continuum Emission 

3.1. Total Intensity 

Figure [T] shows the distribution of the radio continuum emission in NGC 4214. The 
emission is most extended at 20 cm despite the higher noise levels in the 20 cm image. At 
3 cm, the emission is largely confined to two regions near the center of the galaxy. The circles 
in the middle and right panels of Figure Q] show the largest angular scale imaged by the VLA 
at 6 cm and 3 cm. The largest angular scale imaged by the VLA at 20 cm is larger than the 
region shown in Figure [TJ 

We recover most of the radio continuum flux from the galaxy. The first evidence in favor 
of this conclusion is that the extent of NGC 4214 at both 6 cm and 3 cm is much smaller than 
the largest angular scale imaged at those wavelengths. The second line of evidence for this 
conclusion is shown in Figure [2j In this Figure, we plot the sum of the 3 cm radio continuum 
emission in a particular annulus as a function of radius from the center of the galaxy. The 
presence of a large negative bowl surrounding the radio continuum emission would indicate 
that large-scale emission is being resolved out. This effect would be largest for the 3 cm data 
because it will resolve out emission on smaller scales than the 6 cm or 20 cm data. There 
is no significant negative bowl surrounding our 3 cm radio continuum emission. The dip at 
150" is well beyond the main body of radio continuum emission and is the same magnitude 
as the variations at larger radii. Finally, we have plotted the total flux from our data along 
with measurements of the total flux from the literature in Figure [3J We assume that the 
errors on our total fluxes and those of measurements from the literature without quoted 
errors are 20%. This error reflects the fundamental uncertainty in the calibration of radio 
fluxes. Table [3] lists the values and their sources. Our measured radio continuum fluxes are 
larger than the previously measured radio continuum fluxes. Therefore, we are not resolving 
out significant emission on scales larger than our largest angular scale. Our measured fluxes 
are larger than the fluxes from previous observations because all of the previous observations 
were shallow survey observations and these observations are able to detect more faint diffuse 
emission. 

Given the systematic offset between the archival measurements and our data, the het- 
erogeneity of the archival measurements, and the limited number of measurements from our 
observations, we only fit a simple model consisting of a single power law to our radio contin- 
uum spectrum for NGC 4214. We find that the data can be fit by a single power law with a 
spectral indejiB equal to —0.43 ±0.06. This spectral index shows that the spectrum includes 



4 We define the spectral index, a, as S v oc v a . 
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a mix of thermal and non-thermal emission. We see some hint that the spectrum is convex, 
which would suggest that models including free-free absorption of s ynchrotron and /o r free- 
free emission or including synchrotron losses might be appropriate (|Deeg et al.lll993l ). The 
electron densities measured in NGC 4214 by IKobulnicky &: Skillmanl ( 119961 ) are similar to 
the electron densities determined by the modeling of the radio continuu m spectra of othe r 
similar galaxies including free- free absorption of synchrotron emission by lDeeg et al.l ( 119931 ). 
However, our limited data set does not justify the use of more complex models. Additional 
data points at higher and lower frequencies determined using data that matches the spa- 
tial frequency range and sensitivity of our data are necessary to model these more complex 
physical scenarios. 

There are two isolated regions of radio continuum emission located to the north and 
northeast of the main emission region. The first is a region of diffuse 20 cm emission to the 
north of the galaxy at (12 h 15 m 38 s , 36°22'30"), referred to hereafter as the diffuse synchrotron 
region. This region is likely to be part of NGC 4214 based on the presence of the optical and 
Ho; emission in the same location (see §4.11 for further discussion). The bright radio contin- 
uum source with a steep radio continuum spectrum (see § !3.2p at (12 h 15 m 48.9 s , 36°21'54.10") 
is resolved into t wo sources by the Fai nt Images of the Radio Sky at Twenty centimeters 



(FIRST) survey (Becker et al.lll995f ). iMcMahon et al.l (120021 ) is able to identify an opti- 



cal counterpart for this source on the red POSS-I plates; the optical morphology is clas- 
sified as blended (two local maxima with a single set of connected above-threshold pix- 
els). Unfortunately, there are no determinations of the source redshift. The southern- 



most of the two FIRST sources is coincident with x-ray point source 20 from lHartwell et al. 



(120041 ) . The x-ray emission from this source is extremely hard {{H — M)/(H + M) = 0.698, 
(M — S)/(M + S) = 0.231). Based on th e steep spectra index in the radio, the blended source 
identification of IMcMahon et al.l (120021 ). and the hard x-ray colors from the associated x-ray 
source, we conclude that this source is a background AGN. 



3.2. Spectral Indices 



The spectral indices between the 20 cm and the 6 cm data and the 6 cm and the 3 cm 
data are shown in Figure |H Regions in the input images below 3cr were masked prior to 
calculating the spectral index. Thermal emission has a spectral index of -0.1 and synchrotron 
emission has a spectral index of ~ —0.7, although free- free absorption and synchrotron 
losses can modify the observed spectrum. Th e two brightest emissi on peaks are thermally 
dominated with the younger (2.5 to 3.0 Myr; iMacKenty et al.l 120001 ). sout hern peak (NGC 



4214-11) having a flatter spectral index than the older (3.0 to 4.0 Myr; IMacKenty et al. 
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2000l ). northern peak (NGC 4214-1). Flat radio continuum spectra like t hose seen in NGC 



4214-1 and NGC 4214-11 are characteristic of dwarf starburst galaxies (jKlein et al. 1984 ; 
Deeg et al.lll9931 ). The two fainter sources northwest of the galaxy (at 12 h 15 m 34 s , 36°20') 



are also thermally dominated, but with a less steep spectrum than NGC 4214-1 and NGC 
4214-11. The regions between NGC 4214-1 and -II and the peaks to the northwest have steeper 
spectra showing that synchrotron emission plays a more important role in these regions. 



3.3. Polarized Emission 

Since we were interested in investigating the structure of the magnetic field, we deter- 
mined the polarized intensity and the polarization angle from the Stokes Q and U images. 
The der ived polarized intensities were corrected for t he bias in the estimated polarized in- 



tensity (ISimmons fc Stewartlll985l ; IVaillancourtl 120061 ) . 



The percent polarization gives the strength of the uniform field relative to that of the to- 
tal field, while the polarization angle rotated by ~ 90° gives the orientation of the magnetic 
field. The polarized emission of NGC 4214 is shown in Figure [5j There is little polar- 
ized emission associated with NGC 4214. This lack of polarization could be due to several 
instrumental/imaging effects or due to depolarization effects internal to the galaxy. The pos- 
sible instrumental imaging effects are: 1) depolarization due to mosaicing, 2) the polarized 
emission is stronger at lower resolutions, or 3) bandwidth depolarization due to combining 
two IFs. To check that mosaicing the observations has not depolarized the final image or 
produced spurious polarization signatures, we imaged each pointing separately at 6 cm and 
3 cm. The images of each pointing at both frequencies were largely depolarized, so mosaicing 
was not the cause of the depolarization seen in Figure To see what effect decreasing the 
resolution of the images had we made images with 1.5 times and 2.0 times the resolution 
of the VLA images and an image with the resolution of a 100-m single dish telescope (e.g., 
the Effelsberg 100-m telescope or the Green Bank Telescope). Decreasing the resolution 
did not increase the amount of visible polarization, so the resolution of these images was 
not the cause of the depolarization. To see if bandwidth depolarization is the cause of the 
depolarization, we imaged each IF separately and compared their polarized emission and de- 
termined that there was no significant depolarization associated with combining the two IFs 
for imaging. Based on these experiments it does not seem likely that instrumental/imaging 
effects are the cause of the depolarization. 

The lack of polarization in Figure [5] could also be due to the properties of the galaxy 
itself including depolarization effects interna l to the galaxy or ins ufficient sensitivity. Depo- 



larization effects internal to the galaxy, e.g., ISokoloff et al.l ( 119981 ). decrease as ones goes to 
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shorter wavelengths. Therefore, the image at 3 cm should be largely free of depolarization 
effects while still having appreciable synchrotron emission (see Sections 13.11 and 14.3ft . The 
3 cm polarization image does not show appreciably more polarization than the 20 cm polar- 
ization image, so internal depolarization effects do not seem to be leading to depolarization 
at 20 cm and 6 cm. In Section 14. 5. H we use our polarization data to constrain the uniform 
magnetic field strength. 



3.4. Comparison with Previous Observations 



Our data represent the first observations of the polarized radio continuum emission from 
NGC 4214. There are, ho w ever, three previo us stu dies of the radio contin uum emission of 
NGC 4214: kllsoppl jl979h . beck et all (bood ). and lMacKentv et all ( fioooh . 



The earliest interferometric observations of the radio continuum emission from NGC 



4214 were done at 20cm by lAllsoppI (11979I ). Despite their much larger beam (58" by 96") 
and higher noise levels, the image derived from these observations has the same overall 
shape as our higher resolution and better sensitivity image. The emission to the southeast 
at 12 h 15 m 34 s , 36°20' that they suggest is a background object is actually clearly part of the 
main radio continuum distribution in our higher resolution observations and is associated 
with faint Ha and optical emission. 



Beck et all fl2000[ ) observed NGC 4214 with the VLA at 20 cm, 6 cm, 3.6 cm, and 2 cm 
with a resolution of ~ 2". By optimizing their observations for higher resolution , however, 



they resolve out most of the extended structure seen in Figured! iBeck et al.l (120001 ) . however, 
are able to resolve NGC 4214-1 and NGC 4214-11 into several point sources with hints of 
more extended emission, particularly in NGC 4214-1. They find that the spectral indices 
that they derive between 6 cm and 2 cm are rising, i.e., positive, and attribute this to the 
free-free emission becoming optically thick. They estimate that the turnover frequency for 
NGC 4214 is closer to 15 GHz than to 5 GHz. We do not identify any rising spectrum 
sources between 6 and 3 cm. The beam of our lower resolution observations includes both 
optically thick sources and extended emission so we would not expect to be able to identify 
individual optically thick sources. 



The 20 cm and 6 cm radio c ontinuum observa tions in iMacKenty et al.l (120001 ) are lower 
resolution (~ 4") than those in IBeck et al.l (120001 ). but still higher resolution than our ob- 
servations. Again, the higher resolution of these observations compared to ours allows them 
to resolve NGC 4214-1 and NGC 4214-11 into several po int sources with h ints o f extended 
emission, but resolves out much of the diffuse emission. IMacKenty et al.l (120001 ) were able 
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to identify one of the points sources in NGC 4214-1 as a supernova remnant based on its 
spectral index and the properties of the associated optical emission, showing that supernovae 
have had enough time to produce appreciable synchrotron radiation in NGC 4214. 

4. Discussion 

4.1. Comparison of the Radio Continuum Emission with Emission at Other 

Wavelengths 

In Figures [6], [7J and El we compare the distribution of the radio continuum emission 
from NGC 4214 with the distribution of its optical and Ha emission. Figure |6] is a three-color 
image relating the Ha (red), optical U band (blue), and optical V band (green) to the 20 cm 
radio continuum emission. In Figures [7J and El we provide more detailed comparisons of the 
20 cm, 6 cm, and 3 cm emission with the Ha and optical V band, respectively. 

From Figure [71 we can see that the radio continuum at 3 cm closely follows the dis- 
tribution of the Ha emission. This correspondence shows that the 3 cm radio continuum 
emission has a large contribution from thermal emission. As the wavelength of our observa- 
tions increases, we see more diffuse radio continuum emission. In the northern half of the 
galaxy, there is diffuse emission that corresponds to the faint Ha filaments at 12 h 15 m 37 s , 
36°20'45". In the south, there is some diffuse emission that does not closely follow the Ha 
distribution at 12 h 15 m 45 s , 36°18'00" and 12 h 15 m 47 s , 36°18'45"; these regions have a large 
contribution from synchrotron emission. The western tip of the diffuse synchrotron region 
(12 h 15 m 38 s , 36°22 / 30") is associated with Ha emission and some faint 6 and 3 cm emission. 
This part of the region is likely dominated by thermal emission. The remainder of the dif- 
fuse synchrotron region does not have any corresponding Ha emission or any appreciable 
emission at 3 or 6 cm, showing that is it dominated by non-thermal emission. 

Figure El shows that the radio continuum emission is confined to the central regions of 
this galaxy. We see that the radio continuum emission, the Ha emission, and the stellar 
emission are all coincident in a bar-shaped region at the center of the galaxy. The diffuse 
synchrotron region (12 h 15 m 38 s , 36°22'30") has stars associated with it. This correlation leads 
us to believe that this feature is an older star formation region. The noticeably resolved, 
optical object in the center of the region is green in Figure El which implies its emission is 
dominated by optical V band emission and thus belongs to an older stellar population. The 
expected lifetime for the free-free emission is about 10 Myr (approximately the lifetime of 
the most massive stars). The synchrotron emission lifetime is set by the confinement time 
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for cosmic rays, i.e., the time it takes for cosmic rays to diffuse out of a galaxyo For our 
Galaxy, this timescale is on the order of 10 Myr. In Section 14. 6[ we estimate that this time 
scale is also approximately 10 Myr for NGC 4214 (although see the cavaets discussed in 
Section |4TB"1) . There is a delay expected between the synchrotron emission and the free-free 
emission since the first supernovae do not appear until 3.5 Myr after the initial burst of star 
formation. This argument suggests that there can be synchrotron emission in a region after 
the free-free emission has dissipated. Note that in this picture, we are assuming that the 
cosmic rays have not diffused very far from their sources on timescales of 10 Myr. 

The timescales for free-free and synchrotron emission allow us to estimate the age of 
the star formation episode that formed the diffuse synchrotron region. The eastern half of 
the region has formed stars greater than 10 Myr ago so the free-free emission has dissipated, 
but less than 13.5 Myr ago, so the cosmic rays are still confined to the galaxy and producing 
synchrotron emission. The western half has formed stars within 10 Myr and thus still has 
both synchrotron and free-free emission. The lack of emission in diffuse synchrotron region 
at 6 and 3 cm is due to the decrease in synchrotron intensity as the wavelength decreases 
(i.e., as the frequency increases). Assuming a spectral index of -0.7 and using the intensity 
of the diffuse synchrotron emission at 20 cm, the synchrotron emission from this patch would 
be below the 3a detection limit at 6 cm. 



4.2. Fraction of Thermal Emission 



To determine the magnetic field strength, we need to determine how much of the radio 
continuum emission is due to non-thermal (synchrotron) emission and how much is due to 
thermal (bremsstrahlung) emission. Since the electrons that produce the thermal emission 
in the radio are the same electrons that produce Ha in the optical, we use a flux calibrated 
Ha image kindly provided by D. Hunter (IHunter fc Elmegreenll2004j ) to estimate the thermal 
contri bution to the radio emission. We use the approximate expression derived by lHunt et al. 
(120041 ) to convert the Ha fluxes to radio continuum fluxes 



mJy 



1.16 1 + 



n(H+ 



T 



10 4 K 



0.617 
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where f v is the radio continuum flux, n(He + )/n(H + ) is the ratio of ionized helium to ionized 
hydrogen, T is the temperature of the region, v is the frequency of the radio continuum 



5 Sm all distortions in the galactic magnetic field randomly scatter cosmic rays. See 
(|2009| ) for a review of the relevant physics. 
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observations, and fu a is the Ha flux. These expressions are valid for temperatures between 
10,000 K and 20,000 K and densities between 100 cm -3 and 1000 cm -3 . We use a value 



of 0.086 for n(He + )/n(H + ) and a temperature of 10,000 K (IKobulnickv fc Skillmanlll996l ). 
The synchrotron flux is then 

F t h 



nth 



tot 



tot 



(2) 



where F nt h is the flux of the synchrotron (non-thermal) emission, F tot is the total flux ob- 
served, and F t h is the thermal flux calculated using Equation p]). Any values of F t h/F tot 
greater than one or below zero are set to one and zero, respectively. In practice, this affects 
only a few edge pixels. 

We correct the Ha flux for foreground extinction using the value derived from the maps 
of ISchlegel et all (119981 ): Ar ~ Au a ~ 0.058. We do not correct for the internal extinction of 



NGC 4214, which var ies as a function of position in the galaxy (jMaiz-Apellaniz et al.lll998l ). 
Lisenfeld et all (120081 ) find that the internal extinction of NGC 4214 for Ha is between 0.3 
and 0.7 magnitudes. The values of high internal extinction in this galaxy are slightly offset 
from the peaks in the Ha emission. Due to our low resolution (~ 14"), both the peaks of 
thermal radio/Ha emission and the regions of high extinction fall within the same beam. 
Not correcting for the effect of internal extinction could us to underestimate the flux of the 
free-free emission by a factor of 1.3 to 1.9 and thus to overestimate the synchrotron flux 
by factors of 1.25 to 2 for regions that have 0.3 magnitudes of extinction and by factors 
2 to 3 for regions that have 0.7 magnitudes of extinction. The lower limits in both cases 
refer to the flux peaks. Given the value of a used in our calculation of the magnetic field 
(see Section HT31 for details), this will change our magnetic field estimate by 9 to 40%, with 
the mean change closer to 9%. Correcting the free-free emission also has an effect on the 
measured polarization of the emission. If the Ha data is not corrected for extinction, the 
flux of the free-free emission will be an underestimate, so the intensity of the synchrotron 
will be an overestimate. Then the percent polarization of the synchrotron emission, which 
is related to the ratio of the uniform field to the random field, will be an underestimate. 



Our thermal flux estimates as a function of position in NGC 4214 are shown in Figure 



4.3. Magnetic Field Strength 



We can calculate the magnetic field strength using the estimated synchrotron flux from 
Section 14.21 We use the revised equipartition estimate of the magnetic field strength given 
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in Beck fc Krause! (120051 ) 



A7r(l-2a)(K + l)I u E l v +2a (u/2 Cl ) 



l/(3-a) 



(3) 

(—2a — l)c 2 (a)/c4 

where B t is the total equipartition magnetic field in G, a is the spectral index of the emission^], 
Kq is the number density ratio of the protons to the electrons, E p is the rest energy of the 
proton, ci is a constant equal to 6.26428 x 10 18 erg -2 s _1 G -1 , c 2 is a constant tabulated 
on page 232 of iPacholczykl (ll970l ) as c 5 , I is the line of sight through the galaxy, and c 4 
is a function correcting for the inclination of the region with respect to the sky. Since our 
polarization observations shown in Section 13.31 do not reveal a large-scale uniform field, we 
use the version of C4 for the case of a completely random, isotropic field 

« = (I) ' (4) 

In this does not explicitly depend on the inclination, but instead represents an 

average over all possible orientations for the random field. 

We need to carefully consider which frequency has the best estimate of the synchrotron 
flux. In Section 13. 1[ we see hints that the radio continuum spectrum of NGC 4214 may 
be affected by either free-free absorption of synchrotron emission and/or synchrotron losses, 
although we do not have adequate data at higher and lower frequencies to constrain either 
model. Free-free absorption of synchrotron emission would cause the 20 cm synchrotron flux 
to be underestimated. Cosmic rays also lose energy via synchrotron emission if they are not 
reaccelerated. This effects the higher energy cosmic rays more than the lower energy cosmic 
rays and would thus cause the 3 cm synchrotron flux to be underestimated. Therefore, the 
6 cm image represents the best wavelength for us to determine the magnetic field strength of 
NGC 4214. We use a typical value for a of -0.7. Given the uncertainties in the mechanisms 
responsible for shaping the radio continuum spectrum of NGC 4214, it is better to use a 
typical value than a value derived from the data. We note that the spectral index between 
the 20 cm and 3 cm synchrotron fluxes for our data is -0.6, which is very close to the typical 
value. 



We can use the following expression from lOtt et al.l (120011 ) to estimate the scale height 
of NGC 4214's disk from neutral hydrogen observations and obtain the line of sight distance 
through NGC 4214's disk 

-1 



h [pc] = 5.79 x 10 21 ( 
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(5) 



6 We use the opposite sign convention for a as 



Beck fc Krausel(l2005l) . 
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where h is the scale height, a gas is the velocity dispersion of the neutral gas, Nhi is the 
column density of the neutral gas, an d pml Ptat is the ratio of the density of neutral gas to 
the total density. From 



Al 



soppI (Il979f ). we know that M HI /M tot is 0.3, so since M m /M i 



tot 



PmlPtot = 0.3. lAUsoppI (Il979f ) also gives us a column density for the neutral gas of 8 x 



10 cm and a velocity dispersion of 10 kms . These values give a scale height for this 



disk o f 200 pc, which agrees with the estimate of the disk scale height by lMaiz-Apeilaniz et al. 
(119991 ) . Given the uncertain but close to face-on inclination of NGC 4214 we take % = 30° 



The scale height of the magnetic field may be larger than the scale height of the HI. In the 

) , while the scale 



Milky Way, the scale height of the HI is 1.6 kpc ( Marasco fc Fraternah] 201 



height of the synchrotron emission is 4.6 kpc at the solar circle (jFerrierdl200ll ). Unfortunately, 



the face-on orientation of NGC 4214 and the lack of information about its overall field 
topology prevent us from determining the scale height of the synchrotron emission from first 
principles. We can estimate the effect of a much larger magnetic field scale height on the 
magnetic field strength estimates for NGC 4214 using the difference between the Milky Way 
HI and synchrotron scale height and Equation ([3]). If the synchrotron scale height of NGC 
4214 is three times larger than the HI scale height, then the magnetic field strength estimate 
above will be 1.3 times bigger than the actual magnetic field strength. 

Our estimate of the magnetic field strength in NGC 4214 is given in Figure [10J The 
magnetic field strength peaks near the center of the galaxy with a value of 30 pG and 
decreases to a value of 10 pG at the edge of the galaxy. To get an idea of the uncertainties 
in these estimates, we have varied K to 40, a by 0.1, changed the estimated synchrotron 
emission by a factor of 2, changed the line of sight distance through the galaxy by 20%, and 
changed the distance by 20%. With these uncertainities in the initial parameters, we get 
an uncertainty for the central magnetic field strength of 9.5 pG and an uncertainty for the 
magnetic field strength on the edges of 3 pG. 



4.4. Importance of the Magnetic Field in the ISM of NGC 4214 

To determine the relative importance of the magnetic field in the larger context of 
the interstellar medium of NGC 4214, we calculate total magnetic pressure as well as the 
magnetic pressure in the inner and outer regions of NGC 4214 and compared this pressure 
to the thermal pressure from the hot ionized gas, the thermal and turbulent pressure in HII 
regions, and the gravitational "pressure" . These values, along with the values and equations 
used to calculate them, are given in Table HJ 

From these estimates, we can see that the total hot gas pressure, the total magnetic 
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pressure, and the gravitational "pressure" are the same order of magnitude, with the hot gas 
having the largest pressure and the magnetic field the lowest The thermal and turbulent 
pressures of the HII regions, however, dominate the pressure balance i n the central star 
forming regions (although not the galaxy as a whole). IWilcots fc Thurowl ( 120011 ) found that 
the energy in the ionized gas was comparable to that injected via stellar winds. The ISM 
pressure balance of NGC 4214 is similar to that of NGC 1569 JKeplev et all bold ): the 
turbulent and thermal pressures of the HII regions dominate the pressure balance of NGC 
1569 in the central star forming regions, while the magnetic field, hot gas, and gravity are 
approximately equal. Although the magnetic field strength in both galaxies is comparable 
to that of larger spiral galaxies, neither NGC 4214 nor NGC 1569 appear to be magnetically 
overpressured compared to their gravitational potential. The magnetic field of NGC 4214 



may play a role in shaping any outflow of gas from this galaxy as in NGC 1569 (IKepley et al. 



20101 ) . but given the orientation of NGC 4214, we are not able to determine this with our 



observations (see Section 14.5.31) . 



4.5. The Structure of the Magnetic Field of NGC 4214 

4-5.1. Uniform Field Using Polarization Information 



Using our data, we can calculate an upper limit on the uniform magnetic field strength 
in NGC 4214. One can approximate the percent polarization as 



P 
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(6) 



where p is the percent polarization, P is the intensity of the polarized emission, / is total 
intensity of the synchrotron emission, p is the intrinsic polarization of synchrotron emission 
(~ 0.75), is the uniform magnetic field in the plane of the sky, and B t ± is the total 
magnetic field in the plane of the sky ( IBurnl Il966t iBeck et al.l 120031 ). For NGC 4214, we 
calculated the total field in Section 14.31 assuming that the field was entirely random, so 
Equation ([6]) becomes 

Bl , 



P = Po 



2R2 
3 n t 



(7) 



Assuming a value of 40% thermal emission at 6 cm (see Section l4~2l) . we would detect 
at the 3a level 30% polarization at 12 times the noise level (third contour in Figure [5]). In 



7 Since the gravitational pressure was estimated using the Virial Theorem, it should only be compared to 
the pressures for the entire system, not the pressures for individual regions. 
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Section 14.31 we estimate the total magnetic field strength at the 3rd contour to be 15 
Using Equation (J7|), we then get an upper limit on the strength of the uniform component 
of the magnetic field in the plane of the sky of 8 /zG. The large thermal component of the 
radio continuum emission of NGC 4214 reduces our sensitivity to the uniform field. If the 
radio continuum emission in NGC 4214 was all synchrotron, our upper limit on the uniform 
magnetic field strength would be 6 /xG. 

From this analysis, we conclude that we do not detect a uniform component of the 
disk magnetic field of NGC 4214 greater than 8 fiG on size scales larger than our resolution 
(200 pc). The disk of NGC 4214 may have a uniform field, but it would have to be weaker 
than 8 fiG. There could be structure in the magnetic field on size scales smaller than 200 pc, 
but our observations cannot resolve these structures. 



4-5.2. Could the Shape of the Diffuse Synchrotron Region Be Due to a Uniform Field? 

As a star-forming region ages, the relativistic electrons diffuse away from their birth- 
places along magnetic field lines. In a completely random field, one would expect that the 
electrons would diffuse outward isotropically. If there is a uniform component to the field, 
however, the relativistic electrons will preferentially diffuse along the field lines rather than 
across them. The diffuse synchrotron region in NGC 4214 is elongated along a position angle 
of about 60 degrees, which suggests that there is a uniform magnetic field component in that 
direction. 

We can use the distance between the stars and the end of the region (490 pc) and the 
probable age of the region (10 Myr) to estimate the speed at which the relativistic electrons 
diffused. This gives a lower limit on the diffusion speed of 4.7 x 10 6 cms -1 . Assuming that 
relativistic electrons are moving at the Alfven speed (va) 

.2 & 



Airp 



where B is the magnetic field strength, p is the density, and all quantities are in cgs units, if 
we have an estimate of the density in the region then we can get a lower limit for the strength 
of the magnetic field responsible for diffusing the cosmic rays. The speed of the cosmic rays is 
coupled to the density of the thermal electrons, so we can use the Ha image described above 
to get an estimate of the density of the thermal electrons in the diffuse synchrotron region. In 
the diffuse synchrotron region there is 4.11 x 10~ 6 erg cm -2 s _1 sterad -1 of emission. From 



Recall that this field strength was derived in Section 14.31 assuming that the field was entirely random. 
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Spitzeri fll978f ) 



EM [pc cm" 6 ! = $ Iu ^ (9) 

LP J 2.46 x 10 17 /ica 3j2 

where EM is the emission measure, I v is the specific intensity, 0:3.2 is the production coefficient 
for the Ho line, and all values (unless noted) are in cgs. We use a value of a 3j2 of 11.7 x 10 -14 
which is valid for 10,000 K. This gives us an emission measure for the region of 47 pc cm -6 . 
If we assume that NGC 4214 is nearly face-on, we use the line of sight distance through the 
disk derived in Section 14.31 to convert the emission measure to a number density of thermal 
electrons. This calculation gives a number density for the electrons of 0.34 cm -3 . Assuming 
that the ionized gas has the same number of electrons and protons gives a density for the 
ionized gas of 5.7 x 10~ 25 gem -3 . The lower limit on the total magnetic field strength can 
then be derived using 

B = VA^Tp. (10) 

Using our values of va and p, the lower limit on the total magnetic field strength in the diffuse 
synchrotron region is 13 /zG. This field strength is the same as the field strength derived for 
the diffuse region from the 20 cm data using the assumptions described in Section 14.31 At 
20 cm, we are sensitive to polarizations greater than 26% in the diffuse synchrotron region. 
Using the magnetic field strength determined from the synchrotron emission (13 piG), we get 
an upper limit on the uniform field strength of 7.6 p,G. 

We can use the ratio of the length of the region to the width of the region to derive 
the ratio of the uniform magnetic field strength (B u ) to the random magnetic field strength 
(B r ) and compare it to the limits on B u j B r from our data. Let us take y to be the direction 
along the diffuse synchrotron region and x to be the direction across it, the equation of the 
magnetic field to be B = yB y + xB x , X to be half the total width of the region, and Y to be 
half the total length of the region. B y and B x are both constant in a correlation length Z, but 
B x can be positive or negative. In this picture, B y corresponds to the uniform component 
and B x corresponds to the random component. Then the equation of a field line is 

? = ^- (11) 
dx B x 

The displacement in y when the field line has gone a distance I in x is 

Ay = ^. (12) 



'X 



The cosmic rays stream at the Alfven speed, so the time to travel the distance / is r = 1/vax, 
where va x — B x /{A^p) 1 / 2 . After a time t, the cosmic rays have taken N = t/r random 
steps of length I in the x direction , so the rms displacement of the cosmic rays, X, is 
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IN 1 / 2 = {IvAxt) 1 / 2 ■ The displacement in the y direction is just NAy 
elongation of the cloud is 



Y 
X 



VAxt 
I 



1/2 



By 

B„. 



N i/2^y 
B' 



NlB y /B x . Then the 
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We can estimate N from the strength of the random field derived above from our data 
(12.9 /xG), the density of the region (0.34 cm" 3 ), the age of the region (10 Myr), and the outer 
length scale for turbulence (90 pc in the LMC; iGaensler et al.l 120051 ) . This gives a value for 
N of approximately 6. The measured values of X and Y are 235 pc and 491 pc, respectively, 
which gives a value for Y/X of 2.1. With our estimate of N, the ratio B y /B x = B u /B r is 
0.85. From our data, the ratio of B u /B r allowed by our data is 0.69, which is close enough 
to the ratio estimated from the size of the region to make our scenario that the diffuse 
synchrotron region is elongated by a uniform field plausible. 



4-5.3. Possible Structure in the Halo of NGC 4^14 



The analysis in Sections 14.5.11 and 14.5.21 only applies to the structure of the magnetic 
field in the disk of NGC 4214. The nearly face-on orientation of NGC 4214 does not allow 
us to determine whether NGC 4214 has a uniform mag netic field oriented along a galactic 
outflow like that seen in NGC 1569 (IKepley et al.ll2010l ). From estimates of the mechanical 
energy in the wind (L mer h ~ 10 40 — 10 41 er g s _1 ), it appears that it is possible for NGC 4214 



to experience a blow-out (jOtt et al.ll2005bl ). Therefore, NGC 4214 might have magnetic field 
structure perpendicular to its disk in its halo. Detection of this halo field may be possible 
using t he dense grid of backg round rotation measures produced by the Square Kilometer 
Array (IBeck fc Gaenslerl 120041 ) . 



We would like to note that while in principle the rotation measure of the background 
source to the northeast of the main body of emission could be used to estimate the magnetic 
field direction and strength along the line of sight, i.e., through the disk, unfortunately there 
is no Ha emission in the region of the galaxy associated with the background source, so 
we would be unable to determine a field strength. In addition, there are no other nearby 
background sources to allow us to estimate an average rotation measure for background 
sources in that region of the sky. 
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4.6. Cosmic Ray Lifetimes 

Using the arguments from Section 14.5.21 we can also estimate a timescale for cosmic ray 
confinement in NGC 4214. Assuming that a supernova goes off at the midplane of the disk, 
the cosmic rays would have to travel 200 pc to to escape the disk. Since cosmic rays diffuse 
at the Alfven velocity, the total displacement of the cosmic rays, L, in a time t is 

L - ^ - < 14 > 

The timescale for the cosmic rays to escape the disk is then 

t= IB ' (15) 

where n e is the number density of electrons and mn is the mass of a hydrogen atom. We 
estimate n e to be 1 c m" 3 , B to be 20jzG, and take I to be the outer length scale of 



the turbulence, which iGaensler et al.l (120051 ) found to be 90 pc in the LMC. Using these 
parameters, the timescale for the diffusion of cosmic rays out of the disk of NGC 4214 is 
approximately 10 Myr - comparable to the diffusion time in the Milky Way. If the scale 
height of the synchrotron emission is much larger, as discussed at the end of Section I4.3[ 
this timescale cou ld be as long as 67 Myr. NGC 4214 may have a wind i n its central portion 



Ott et al.l (l2005bl ) which would shorten the cosmic ray confinement time (jEverett et al.ll2008 



20101 ). Detailed modeling of the effect of a central wind on the cosmic ray confinement time 



in NGC 4214 is beyond the scope of the paper. 

Synchrotron losses will further limit the lifetime of synchrotron emission. To see if the 
synchrotron loss timescale will dominate over the cosmic ray confinement time, we calculate 
the lifetime of cosmic ray electrons due to synchrotron losses. The power P in erg s -1 
radiated by synchrotron emission is 

P = ^rcPVf (16) 

where ot is the Thomson cross-section, c is the speed of light, j3 is v/c, and 7 is (1 — (3 2 )~ l l 2 . 
Since cosmic rays are traveling very close to the speed of light, we can assume /3 ~ 1. 
Assuming that all the emission is produced at the critical frequency 

v c = —7 — sin a (17) 
Air mc 

and using the relation E = ^mc 2 , we find that the energy of the electron at the critical 
frequency is 

/ 47rm^c 5 z/ r .\ ^ 2 
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where the factor of sin a has been ignored since it is a factor on the order of one. Dividing the 
result of Equation f fl8|) by Equation f fl6|) gives an estimate for the synchrotron loss timescale. 
For emission at 20 cm (1.4 GHz) and a magnetic field strength of 10 /zG (the field strength 
in the diffuse synchrotron region), the synchrotron loss timescale for a cosmic ray electron 
is 44 Myr. At 6 cm (8.5 GHz), the lifetime of the cosmic rays due to synchrotron losses is 
shorter - 18 Myr. Therefore, at 20 cm in the diffuse synchrotron region, the confinement 
time is shorter than the synchrotron loss timescale and thus the confinement time dominates. 
(Although if the synchrotron scale height is 600 pc rather than 200 pc, the synchrotron loss 
time will dominate.) At lower frequencies and larger magnetic fields, the synchrotron loss 
time scale decreases so that in the center of the galaxy with a magnetic field of 30 fiG and 
at 3 cm, the synchrotron loss timescale may be as short at 6 Myr. 



4.7. Generating the Magnetic Field 

To evaluate the possibility of an a — u dynamo operating in NGC 4214, we compare 
the shear and turbulence in NGC 4214 to that of a larger galaxy, where ana-w dynamo 
may be operating: the Milky Way. 

First, we evaluate the effect of differential rotation by calculating the shear in both 
systems. Shear is defined as dQ/dr, which for a flat rotation curve is just v/r 2 . At the solar 
radius (8 kpc), the rotation speed is 220 kms -1 , so the shear in the Milky Way is about 
3.4 kms~ 1 kpc -2 . In Figure [TT| we show the absolute value of the shear in NGC 4214 as a 
function of radius. To calculate th e shear, we used the neutral hydrogen rotation curve for 



NGC 4214 given in lAllsoppI ( 119791 ). corrected for the distance and inclination used in this 
paper. The shear in NGC 4214 varies significantly below a radius of 2 kpc, but the average 
value at these radii is close to the Milky Way value. At a radius of 2 kpc, the shear begins 
to increase, reaching a maximum of 12 kms _1 kpc~ 2 at a radius of about 3 kpc. The shear 
then begins to decrease and by 5 kpc is at or below the shear estimated at the solar radius in 
the Milky Way. We would note that the diffuse synchrotron region discussed in Section l4.5.2l 
is located in the region of maximum shear about 2.5 kpc from the center of NGC 4214. 
According to Rayleigh's Criterion, the quantity Qr 2 must increase as a function of r for 
a rotation curve to be stable. This is satisfied for the Q shown in Figure [11] everywhere 
except in 3 kpc < r < 5.5 kpc, but the observational uncertainties are large enough that the 
departure from stability, which is slight, is most likely due to error. 

To evaluate the turbulent component of the dynamo process, we assume that the turbu- 
lence in these galaxies is generated by supernovae and use the supernova rate as a proxy for 
turbulence. In the Milky Way, the supernova rate from type II supernovae is 27 Myr -1 kpc~ 2 
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rate. 



(IFerrierell200l[). For NGC 4 2 14, we can determine the supernova rate from the star formation 



Hunter fc Elmegreenl (120041 ) determine a star formation rate per unit area for NGC 
4214 of 5.0 x 10 3 M Q yr~ 1 kpc -2 . Assuming a Salpeter IMF, that the range of masses for 
stars is 0.8 M to 150 M , and that only stars above 8 M become Type II supernovae, we 
obtain a rate for the formation of stars greater than 8 M of 85 Myr -1 kpc~ 2 . Assuming 
a massive s tar lifetime of 30 Myr and that the burst of star formation in NGC 4214 began 
4 Myr ago (IMacKenty et al.ll2000l ). that gives a supernova rate of 11.3 Myr -1 kpc -2 . 



Comparing the shear and the supernova rates between the Milky Way and NGC 4214, 
we see that NGC 4214 has a similar level of shear, but a lower supernova rate by a factor of 
2. This rough calculation may indicate that additional turbulence needs to be injected by 
supernovae in order to provide the turbulence necessary for the shear to generate a large- 
scale field. Given that the episode of star formation in NGC 4214 started relatively recently 
and that it seems to be continuously producing stars, we anticipate that supernovae will be 
able to generate additional turbulence as the galaxy ages. In shear-dominated turbulence, 
one expects that the field would be oriented (although not necessarily directed) azimuthally. 
The region of diffuse synchrotron emission is located near the shear maximum in NGC 
4214, which suggests that the magnetic field in this region is oriented by the shear, but not 
necessarily directed as it would be in a large scale dynamo. 



Amplification by an a — u mechanism is rather slow, so various authors (e.g. iParker 



19921 ) have proposed mechanisms to decrease the time needed to amplify the field. If a large- 
scale field c ould be generated for NGC 4214 u sing a faster dynamo like the one proposed for 
NGC 4449 (jOtmianowska-Mazur et al.ll2000l ). it may have not had time to develop yet or 
it may have been disrupted by recent star formation. Note that we cannot probe magnetic 
field structure on size scales less than 200 pc due to the resolution of our observations. 



5. Conclusions 

In this paper, we have presented the deepest radio continuum polarization images of 
the irregular galaxy NGC 4214 to date. Below we summarize the overall conclusions of our 
work. 

The global radio continuum spectrum for NGC 4214 is convex, which suggests that either 
free-free absorption of synchrotron emission and/or free- free emission or synchrotron losses 
are occurring. We fit the global radio continuum spectrum of NGC 4214 with a single power 
law with a spectral index of -0.43 because of the offsets between our data and the archival 
data, the limited wavelength coverage of our data, and the heterogenity of the archival data. 
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This spectral index indicates that the radio continuum emission in NGC 4214 is a mix of 
thermal and non-thermal emission. The spectral indices of the radio continuum peaks are 
flat, with the younger region (NGC 4214-11) having the flatter spectral index. Between the 
radio continuum peaks, the spectral indices are steeper, showing that synchrotron emission 
plays a more important role there. We do not detect any significant polarized emission in 
NGC 4214 at size scales greater than 200 pc and conclude that this is due to the intrinsic 
nature of the galaxy rather than instrumental/imaging effects or internal depolarization. 

We find that the radio continuum emission at 3 cm is correlated with the Ha emission 
from this galaxy, while the radio continuum emission at 20 cm is more extended. The central 
regions of NGC 4214 are approximately 50% thermal at 3 cm. North of the main body of 
radio continuum emission, we detect a region of diffuse synchrotron radiation associated 
with a group of stars in NGC 4214 and speculate that this region is a relic of previous star 
formation. Based the relative timescales of free-free and synchrotron emission, we estimate 
an age for this region of between 10 and 13.5 Myr. 

We estimate the fraction of thermal radio continuum emission in NGC 4214 using Ha 
images. Using this estimate, we calculate the estimated intensity of the synchrotron radiation 
to estimate the magnetic field strength of NGC 4214. The field strength reaches a maximum 
of 30±9.5 /xG in the center of the galaxy and tapers off to 10±3 /iG at the edges. Comparing 
the pressures of the various components of the ISM of NGC 4214, we find that the hot 
gas pressure, the magnetic pressure, and the gravitational "pressure" all the same order of 
magnitude. In the central star forming regions of NGC 4214, the thermal and turbulent 
pressure of the HII regions dominate the pressure balance of this galaxy. The pressure 
balance in NGC 4214 is similar to that of another irregular galaxy: NGC 1569. In both 
galaxies, star formation plays the dominant role in shaping the ISM. Although the magnetic 
field strength in both galaxies is comparable to that of larger spiral galaxies, neither NGC 
4214 nor NGC 1569 appear to be magnetically overpressured compared to their gravitational 
potential. 

Using our polarization data, we place an upper limit on the uniform magnetic field 
strength of 8 /xG on size scales greater than 200 pc. We speculate that the elongation of the 
diffuse synchrotron emission is due to a uniform field component. Using simple estimates of 
the Alfven speed and the density of the region along with upper limits on the polarization 
from our data, we place an upper limit of 7.6/xG on the uniform field strength in this region 
and calculate a total field strength of 13 /xG and a random field strength of 11 fiG. These 
limits are consistent with the observed elongation of the diffuse synchrotron region. 

We use the simple model developed in Section 14.5.21 to estimate the cosmic ray confine- 
ment time in NGC 4214 and find it to be appropriately 10 Myr, comparable to the diffusion 
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time in the Milky Way. Synchrotron losses and a galactic wind may further limit the lifetimes 
of cosmic rays in the galaxy. 

Finally, we investigate whether an a — u dynamo could operate in NGC 4214 by com- 
paring its shear and supernova rate to that of the Milky Way. We find that the shear in 
NGC 4214 is comparable to the Milky Way, but that the supernova rate is half that of the 
Milky Way. We suggest that there is not enough turbulence yet in the NGC 4214 to drive 
ana-w dynamo, but that as the star formation in NGC 4214 progresses the additional 
turbulence will allow ana-w dynamo to operate. 

Of the low mass irregular galaxies with measured magnetic fields, we find that the 
magnetic field of NGC 4214 most resembles that of IC 10, i.e., strong with some evidence of 
uniform fields. In a future paper, we will investigate the properties of all observed magnetic 
fields in irregular galaxies and attempt to create an integrated picture of how they are 
generated and sustained in these galaxies. 
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A. Deriving the Gravitational Pressure 

To derive the gravitational pressure, we start with the Virial Theorem 

2K + U = (Al) 

where K is the kinetic energy and U is the potential energy. Here we assume that the 
density distribution of NGC 4214 can be approximated by an isothermal sphere with density 
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distribution 



p(r) = Po 



r 



(A2) 



(IKepley et al.l 120071 : Ivan Eymeren et al.l 120091 1 . where po and tq are the central density and 
core radius, p is density, and r is radius. Note that the rotation curves of galaxies are more 
typically fit with a pseudo-isothermal sphere 

-i -l 



P(r) = Po 



1 + 



roV 
r J 



(A3) 



which does not have infinite density at zero radius like Equation ( 1A2I) does. As r becomes 
large, the density from Equation flA2[) tends to the density from Equation ( 1A3j) . We are 
measuring the potential far from the center of the galaxy, so we assume that Equation ( 1A2j) 
applies, which yields a much more analytically tractable result. The mass interior to radius 
r, M(r), is then just 



M(r) = / Ati p(r')(r') 2 &r' 
Jo 

The potential energy of an individual shell, dU, is then 

-GM{r) dm 



dU 



-GM(r) Anr 2 p{r) dr 
r 

-167i 2 Gp 2 r*r 



(A4) 
(A5) 

(A6) 

(AT) 
(A8) 



where G is the gravitational constant. Integrating this expression out to a radius R and 
substituting for p gives the potential energy, U, for an isothermal sphere 

~ GhP (A9) 



u 



R 



For an ideal gas, the kinetic energy is just 



K 



-PV 



(A10) 



where P is the pressure and V is the volume of the galaxy. Inserting the expressions for 
the potential and kinetic energies into the Virial Theorem and dividing by the Boltzmann 
constant (k) to put the pressure in units of Kcm -3 , we obtain the expression 



P 
J 



GM 2 

mvk' 



(All) 
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Since V = (4/3)nR 3 and the circular velocity v at radius R in the galaxy is just v 2 = GM/R, 
Equation (lAlip can be rewritten as 



We would note that regardless of the potential equation we use (isothermal sphere, 
constant density, p oc r 2 , etc.), the multiplicative factor in front of the potential changes 
by only a factor of 2. Since the gravitational pressure is proportional to the total potential 
energy, changes in the assumed form of the potential only change the gravitational pressure 
by a factor of a few. 
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Table 1. Summary of VLA Observations 



Observing Band 20cm 6cm 3cm 



Array 


C 


D 


D 


Date 


2005 Aug 7 


2005 Dec 5 & 19 


2005 Dec 


IF 1 Frequency (GHz) 


1.3649 


4.8851 


8.4351 


IF 2 Frequency (GHz) 


1.4351 


4.8351 


8.4851 


Bandwidth/IF (MHz) 


50 


50 


50 


Number of Channels/IF 


1 


1 


1 


Primary Calibrator 


1331+305 


1331+305 


1331+305 


Secondary Calibrator 


1227+365 


1146+399 


1146+399 


Number of Pointings 


1 


4 


16 


Field of View (') 


30 


9 


5.4 


Resolution (") 


12.5 


14.0 


8.4 


Largest Angular Scale (') 


15 


5 


3 


Integration Time per Pointing (hours) 


10.0 


3.73 


1.2 



Table 2. Final Images 



Frequency 


Beam 


PA 


°7 


CT Q 


<ru 


GHz 


// 


O 


fijy beam -1 


fijy beam -1 


fijy beam -1 


1.397 


14.18" x 11.46" 


-77.3° 


48.8 


13.0 


12.9 


4.860 


14.18" x 11.46" 


-77.3° 


25.1 


17.8 


18.5 


8.460 


14.18" x 11.46" 


-77.3° 


24.1 


20.7 


21.0 
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Table 3. Radio Continuum Flux Spectrum 



frequency 


■ni 

Flux 


Mux brror 






GHz 


mJv 


mJv 


Telescope 


Source 


1.40 


51.5 


10.3 


VLA 


6 


1.40 


38.3 


7.7 


VLA/NVSS 


5 


2.38 


36.0 


3.0 


Arecibo 


4 


4.86 


34.0 


6.8 


VLA 


6 


4.85 


30.0 


4.5 


GB 300ft 


2 


4.85 


30.0 


7.0 


GB 300ft 


3 


8.46 


20.5 


0.5 


VLA 


1 


8.46 


24.2 


4.8 


VLA 


6 



1991 



1978 



Referen ces. - (1) ISchmitt et al.ll2006c (2) iBecker et al. 



f3) iGreeorv fc Condon 



(5) ICondon et alJbooi (6) This 



1991 



(4) 



Dressel fc Condon 



paper 



Table 4. Pressures of Various Components of the ISM of NGC 4214 



Pressure 



Component 


Input Values 


Equation 


Reference 


10 5 K cm" 3 


Magnetic field (center) 


B = 30 fiG 


B 2 /(87rfc) 


1 


2.6 


Magnetic field (edges) 


B= IO/jG 


B 2 /(87rfc) 


1 


0.3 


Magnetic field (total) 


B = 13/xG 


B 2 /(8nk) 


1 


0.5 


Hot gas (central) 


n e = 0.147cm" 3 , T = 2.80 x 10 6 K 


2n e T 


2 


8.2 


Hot gas (edges) 


n e = 0.033 cm" 3 , T = 2.68 x 10 6 K 


2n e T 


2 


1.7 


Hot gas (total) 


n e = 0.064cm" 3 , T = 2.17 X 10 6 K 


2n e T 


2 


2.8 


HII regions (thermal) 


n e = 100 cm" 3 , T = 10, 500 K 


2n e T 


3 


10.5 


HII regions (turbulent) 


n e = 100 cm" 3 , v tur f, = 50kms —1 


n e m H v 2 tur J{k) 


1 


300 


Gravity 51 


v rad = 20kms" 1 ,i = 30°, r = l.llkpc 


(v rad /smi) 4 /(4-Kr 2 Gk) 


5 


1.9 



a See Appendix fX] for a derivation of the expression used to calculate the graviational pressure. 



References — (1) This work; f2t lOtt et all ll2005bh ; (3) iKobulnickv fc Skillmanl <1996h ; (4l lWilcots fc Thurowl J200ll^ ; (5) 
lAllsoprj h979h . 
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Fig. 1.— Maps of the total intensity radio continuum emission of NGC 4214 at 20 cm 
(left panel), 6cm (middle panel), and 3cm (right panel). The images are all on the same 
logarithmic scale stretching from 72.3 /zJy beam -1 to 20 mJy beam -1 . The contours are 3, 
6, 12, 24, 48, 96, 192, and 384 times the la noise level given in Table |2j The beam is boxed 
in the lower left corner of each panel and is identical for all three images. The diameter of 
the circles in the middle and right panels show the largest angular scale imaged by the VLA 
at 6 cm and 3 cm, respectively. The largest angular scale imaged by the VLA at 20 cm is 
larger than the size of the region shown. 
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Fig. 2. — 3 cm radio continuum emission per 14" annulus about the center of the galaxy as 
a function of radius. The presence of a large negative bowl of emission, which is not seen 
here, would indicate that significant flux is resolved out. 
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Fig. 3. — Radio continuum spectrum of NGC 4214. The open points are previous measure- 
ments and the solid points are our measurements from this paper. The line is the best power 
law fit to our measurements; see the discussion in § 13.11 
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Fig. 4. — Spectral index between the 20 cm and 6 cm data (left panel) and between the 
6cm and 3cm data (right panel). Pixels below 3a in the input images were masked before 
calculation of the spectral indices. The contours show the total intensity radio continuum 
emission at 6 cm. The contours are 3, 6, 12, 24, 48, 96, 192, and 384 times the la noise level 
given in Table [2j The beam is boxed in the lower left corner of each panel. 
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Fig. 5. — Polarized intensity vectors (E-vectors) overlaid on contours showing the total 
intensity radio continuum emission of NGC 4214 at 20 cm (left panel), 6 cm (middle panel), 
and 3 cm (right panel). The polarized intensities have been bias corrected and only vectors 
with a signal to noise greater than 3 are shown. The contours are 3, 6, 12, 24, 48, 96, 192, 
and 384 times the la noise level given in Table [2j The beam and the polarization vector 
scale bar are boxed in the right hand corner of each panel. The beam is identical for all three 
images. The length of the polarization vector scale bar corresponds to a polarized intensity 
of 5.5 x 10~ 5 Jy beam -1 . 
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Fig. 6. — 20 cm total intensity radio continuum emission contours overlaid on a three color 
optical image of NGC 4214. The optical images are courtesy D. Hunter and red represents 
Ha, blue the U band optical, and green V band optical. The contours are 3, 6, 12, 24, 48, 
96, 192, and 384 times the la noise level given in Table El 
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Fig. 7. — Contours of total intensity radio continuum emission at 20 cm (left panel), 6 cm 
(middle panel), an d 3 cm (right panel) overlaid on an Ha image of NGC 4214 kindly provided 
by Deidre Hunter (jHunter fc Elmegreenl 12004 ). The contours are 3, 6, 12, 24, 48, 96, 192, 
and 384 times the la noise level given in Table El The beam is boxed in the lower left corner 
of each panel. 
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Fig. 8. — Contours of total intensity radio continuum emission at 20 cm (left panel), 6 cm 
(middle panel), and 3cm (right pa nel) overlaid on an optical V band image of NGC 4214 
kindly provided by Deidre Hunter (jHunter &: Elmegreenll2006l ). The contours are 3, 6, 12, 
24, 48, 96, 192, and 384 times the la noise level given in Table [2j The beam is boxed in the 
lower left corner of each panel. 
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Fig. 9. — Estimated thermal fraction of NGC 4214 with contours of total intensity radio 
continuum emission at 20cm (left panel), 6cm (middle panel), and 3cm (right panel). The 
contours are 3, 6, 12, 24, 48, 96, 192, and 384 times the la noise level given in Table [2j The 
beam is boxed in the lower left corner of each panel. 
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Fig. 10.— Magnetic field strength (in fiG) calculated using the estimated synchrotron 
emission at 6 cm. The magnetic field has not been calculated for the background AGN at 
(12 h 15 m 48.9 s , 36°21'54.10") since it does not belong to NGC 4214. The beam is boxed in the 
lower left. 
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Fig. 11. — Absolute value of the shear calculated using the neutral hydrogen rotation curve 
from lAllsoppI (119791 ) . corrected for the distance adopted in this paper and for the inclination 
of NGC 4214, as a function of radius. 



